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ABSTRACT 



In radio astronomy the polarimetric properties of radiation are often modified 



, during propagation and reception. Effects such as Faraday rotation, receiver cross-talk, 

^ ' and differential amplification act to change the state of polarized radiation. A general 



description of such transformations is useful for the investigation of these effects and 
for the interpretation and calibration of polarimetric observations. Such a description 
is provided by the Lorentz group, which is intimately related to the transformation 
properties of polarized radiation. In this paper the transformations that commonly 
■ arise in radio astronomy are analyzed in the context of this group. This analysis is 

then used to construct a model for the propagation and reception of radio waves. The 
implications of this model for radio astronomical polarimetry are discussed. 

On ■ Subject headings: polarization — techniques: polarimetric 



Introduction 



o 

' In radio astronomy transformations occur during the propagation and reception 

■ of radio waves that act to change the state of polarized radiation. Some of these 

transformations, such as Faraday rotation in the ionosphere or the interstellar medium, 
^ . arise from propagation effects that may themselves be of astrophysical interest. 

Others originate from instrumental effects such as differential amplification or receiver 
cross-talk, and have an adverse effect on polarimetric observations. Realistically many 
such effects may be present, each having its own time and frequency dependence, 
and collectively acting to distort measurements of the polarized radiation. The 
interpretation and calibration of these observations may be quite complex, and it is 
useful to have a general context in which to describe these transformations. 

Linear transformations of fully polarized radiation were first investigated by Jones 
(1941), who represented transformations of the two-component transverse electric 
field in terms of 2x2 complex matrices now called Jones matrices. This analysis was 
extended to partially polarized radiation by Parrent &: Roman (1960), who used 
Jones matrices to describe the transformation properties of the coherency matrix. 
Alternatively, both fully and partially polarized radiation may be described by the 
Stokes parameters, and their linear transformations may be represented in terms of 
4x4 real matrices called Mueller matrices (Mueller 1948). 
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In fact these transformations are intimately related to the Lorentz group. This 
relationship arises from the fact that for a plane propagating wave Maxwell's 
equations admit two independent solutions, representing the two orthogonal senses 
of polarization. Thus polarized radiation constitutes a two state system, and the 
linear transformations of all such systems are described by the Lorentz group. This 
relationship has long been known in optics, and different aspects have been discussed 
by many authors (e.g. Barakat 1963, Whitney 1971, Cloude 1986, Pellat-Finet k 
Bausset 1992). Recent work has focused on representations of the Jones and Mueller 
matrices (Opatrny &: Pcrina 1993, Brown & Bak 1995, Han, Kim & Noz 1997). Up 
to a multiplicative constant the set of Jones matrices constitute the group SL(2,C), 
which forms the spin 1/2 representation of the Lorentz group. The corresponding set 
of Mueller matrices constitutes the group S0(3,l) and forms the spin 1 representation 
of the Lorentz group. Finally, the relationship between Jones and Mueller matrices is 
represented through the 2-1 mapping between these two groups. S0(3,l). 

The formulation of the transformation properties of polarized radiation in terms 
of Lorentz transformations affords considerable insight, including the interpretation 
of the Stokes parameters as a Lorentz 4- vector, the classification of transformations 
of this 4-vector as rotations or boosts, and the existence of a polarimetric analogue 
to the invariant interval. In this paper these concepts are reviewed in the context of 
radio astronomical polarimetry. Since both linear and circular bases are widely used 
in radio astronomy, a basis independent formulation is emphasized. This formalism is 
then used to construct a model for the propagation and reception of radio waves. 

2. Representations of Polarized Radiation 

Consider the representation of a transverse electromagnetic wave. Such a wave 
may be described through its two-component transverse electric field vector £{t). This 
vector may also represent the electric field in a waveguide, or the voltage in a pair of 
cables. The vector £(t) is commonly written in terms of the two-component complex 
analytic signal E(t) as £{t) = Re [E(t) expjiwt}]. This construction is familiar from 
both optics (Born & Wolf 1980) and signal processing (Braccwcll 1986). 

For fully polarized radiation the analytic signal is independent of time. In this case 
the relative amplitudes and phases of the two components of E specify the state of 
elliptical polarization of the plane wave. For partially polarized radiation E(t) is time 
dependent, and measurable properties of the wave may instead be defined through 
time-averaging. Such an averaging procedure is conveniently treated through the 
coherency matrix (Wiener 1930, Wolf 1959). This is a 2x2 Hcrmitian matrix formed 
from the direct product of the analytic signal, and may be written as p = (E(t)(8)E'f(i)). 
Here the angular brackets denote time-averaging. As with any Hermitian matrix, the 
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coherency matrix may be written in terms of 4 real quantities {So, S) as 

p = {SoCTo + S • a)/2 (1) 

where (Tq is the 2x2 identity matrix and cr is a 3- vector whose components are the PauH 
spin matrices. These three 2x2 matrices are traceless and Hermitian with determinant 
— 1. The 4 parameters (5o,S) are simply the mean Stokes parameters of the plane 
wave (Fano 1954), with So representing the total intensity. 

Let us now introduce a particular basis. The electric field vector may be represented 
in the Cartesian basis {x,y,z), in which £{t) = {£x{t),£y{t)) is resolved into mutually 
orthogonal components, each orthogonal to the direction of propagation z of the plane 
wave. The corresponding analytic signal is E(t) = {Ex{t), Ey(t)). For the 3-dimensional 
space of S the Cartesian basis {q, ii, v) is used, along with the customary representation 
of the Pauli matrices 



^f^= . r. ] ^0= : ] (2) 






To associate the coherency matrix with the Stokes parameters, we write p in this basis 
as 

From equation |l|, the Stokes parameters in this basis become 

So = {E*it)E^it)) + {E;{t)Ey{t)) Su = 2 Ke[{El{t)Ey{t))] 

S, = {E*{t)EAt)) - {E*y{t)Ey{t)) 5„ = 2 lm[{E*{t)Eym ^ ' 

This is simply the usual definition of the Stokes parameters in a linear basis (Born & 
Wolf 1980). 



3. Transformation Properties of Polarized Radiation 

Let us now consider the transformation properties of polarized radiation. Attention 
is restricted to linear, invertible transformations. This excludes the class of projective 
transformations, which are important in representing perfect polarizing filters. 
Similarly, such transformations cannot describe multipath propagation of coherent 
radiation, such as the focusing or defocusing of radiation by lenses or mirrors. Despite 
this restriction, the set of linear, invertible transformations encompasses a broad 
class of physical processes, including single-particle scattering, propagation through 
anisotropic media, and many transformations arising from instrumental devices. This 
set may also describe the linear transformations of the voltage signal in two cables, 
which are known in linear network theory as two-port networks (Ruston & Bordogna 
1966). This equivalence is particularly useful in radio astronomy, where the two 
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components of the electric field are converted to voltages by a receiver and then passed 
through an electronics downconversion chain. 

The most general linear transformation of the analytic signal may be written as 
E'(t) = tE(t), where the Jones matrix t is a 2x2 complex matrix. As the direct 
product of the analytic signal, the coherency matrix must transform as p' = tp t^. For 
invertible transformations a Jones matrix may be written as t = Vdet t t^v, where 
det t is the determinant of t and t^v is a matrix with unit determinant. The set of 
2x2 complex invertible matrices with unit determinant forms the group SL(2,C), which 
constitutes the spin 1/2 representation of the Lorentz group. 

To investigate the transformation properties of the Stokes parameters, note that 
the determinant of equation ^ is simply det p = — |5p = 5inv (Barakat 1963). This 
is just the form of the Lorentz invariant. Under transformation by the Jones matrix 
t, det p' = l-v/detTpSinv, so that this interval is preserved up to a multiplicative 
constant. The set of transformations that preserve this interval forms the group 
S0(3,l), which constitutes the spin 1 representation of the Lorentz group. That is, the 
Stokes parameters transform as a Lorentz 4-vector, with the total intensity acting as 
the timelike component and the remaining Stokes parameters acting as the spacelike 
components. The condition that the total intensity > restricts this 4-vector to 
lie within or on the surface of the forward light cone. These two cases correspond to 
partially polarized (S'inv > 0) or fully polarized (Sinv = 0) radiation, respectively. 

The representations of the groups SL(2,C) and S0(3,l) are well known in physics, 
but are not widely used in astronomy. Basis-independent representations of these 
groups are now reviewed, and are interpreted in the context of polarimetry. This will 
serve both as an introduction and to establish the notation used in the next section. 
For similar reviews, see Brown & Bak (1995) or Tung (1996). 

The group SL(2,C) contains as a subgroup the set of 2x2 unitary transformations 
SU(2). Any such unitary transformation may be parameterized as 

rn((/)) = e(^''-"'^) =fTo cos0 + i(j-fisin(/> (5) 

where fi is a unit 3- vector. This is called the axis-angle parameterization of SU(2). 
The angle (/> differs from the definition customary in classical and quantum mechanics 
by a factor of 1/2, but is in agreement with the conventions of optics. Under the 
transformation of equation ^, p — > rf^{(p)p rft(— (^) = {aoSo + cr ■ S')/2, where 

S' = Scos20 + S X nsin20-F (n- S)n(l - cos2(/>) (6) 

and we have used the relationship (a ■ a){a ■ h) = a ■ h + ia ■ (ax b). But S' is simply 
the vector resulting from a rotation of S about an axis n by an angle 24> [cf. Goldstein 
1980). This reflects the well-known mapping between SU(2) and the group S0(3), 
whose elements form a representation of the rotations of a 3-dimensional vector. The 
mapping is 2-1, since the two rotations rn((/)) and rn((^ -|- vr) = — r-|-n(0) result in the 
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same vector S'. Such rotations preserve the degree of polarization |S|/5o of the plane 
wave, and are readily interpreted geometrically in the space of the Poincare sphere in 
terms of the axis n and angle 2(f) of rotation. Equation ^ may be represented in the 
basis of equation ^ as 

/ cos (/) + in,, sin (/) + n^,) sin \ ._. 

y [iriu — n^) sva. (f) cos (p — lUg sm (p J 

where n = (nq,n^,n^). From equation ^ the corresponding rotation of the Stokes 
parameters {So, Sg, S„, Sy) in this basis is 

/ 1 



Rfi(2</.) 



rig + {1 — rig) cos 2(f> rigTiu sin^ (j) + sin 2(f> rigriy sin^ (j) — riu sin 2(j) 
n^n^ sin^ (/> — n„ sin 2(/) + (1 — n^) cos 2(/) n^n^ sin^ (/> + sin 2(/) 

yo n^n^ sin^ (/> + n^j sin 2(/) ntjn^ sin^ — sin 2(/) + (1 — n^) cos 20 

(8) 

For example, for {n^,n{^,n{,) = (1,0,0) equation ^constitutes the rotation about the 
q axis Rq(2(/)), which may be interpreted from equation ^ as generating a phase delay 
between the two components of the electric field. 

Next let us consider the group SL(2,C), which has a parameterization similar to its 
subgroup SU(2). Any element of the group may be written as exp (ia • fi + o" • m/3), 
where n and rh are unit vectors. In analogy with equation ^, let us consider 
transformations for which = 0. Such transformations may be written in terms of the 
Hermitian matrix 

t»m(/3) = exp (cr • m /?) = (Tq cosh /3 + a ■ rh sinh (5 (9) 

Under this transformation the coherency matrix becomes p ^rh{P)p ^miP) = 
{(JoS'o + (7 ■S')/2, where 

S'^ = So cosh 2/3 + S • rh sinh 2/3 
S' = S+(5o sinh 2/3 + 2S-msinh2/3)m ^^^^ 

This is simply the result of performing a Lorentz boost on the 4-vector S = {So, S) 
along the axis rh by a velocity parameter 2(3 {cf. Jackson 1975). As in the general 
case, such a transformation preserves the invariant interval S'inv Unlike rotations, it 
does not preserve the degree of polarization |S|/5o of the plane wave. For example, 
there exists some transformation that will completely depolarize partially polarized 
radiation. This is just the analogue of the statement in special relativity that there 
always exists a reference frame in which two events separated by a timelike interval 
occur at the same position in space. Equation ^ may be represented in the basis of 
equation ^ as 

^ _ I cosh /3 + mq sinh /3 (m„ - im^) sinh/3 \ ^^^^ 
\ {rriu + im^) sinh (3 cosh /3 — sinh /? J 
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where rh = {mq,mu,m{,). The boost of equation ^ in this basis becomes 



Ba(2/3) 



/ cosh 2/3 
mq sinh 2/3 
niu sinh 2/3 
y niv sinh 2/3 



mq sinh 2(3 
1 + 2m^ sinh^ /3 
Iniqiriu sinh^ /? 
2mqm^ sinh^ /? 



m„ sinh 2/3 sinh 2/3 ^ 

Imqiriu sinh^ /3 2mqmy sinh^ /3 

1 + 2m^ sinh^ /3 2mum^ sinh^ /3 

2mumy sinh^ /? 1 + 2m^ sinh^ /3 y 



(12) 



The rotations of equations |^ and |8| and boosts of equation 11 and 12 constitute a subset 
of the Lorentz transformations that suffices for the radio astronomical apphcations 
presented in the next section. 

Finally, systems consisting of multiple physical processes are readily modelled 
through successive application of the above transformations. Of particular use in the 
analysis of such composite systems are the commutation relations 



[Rn(a)Rn(/3)] = [Rn(a)BA(/3)] = [Bi,{a)BM] = 



(13) 



As we shall see in the next section, these relationships are useful in determining 
whether the transformations that arise from different physical processes or instrumental 
elements commute with one another. 



4. The Propagation and Reception of Radio Waves 

Let us examine some practical situations in which linear, invertible transformations 
of polarized radiation arise. One such example is a rotation of the electric field vector 
about the direction of propagation. This transformation may represent a rotation of a 
physical device with respect to the plane wave. Another example arises from Faraday 
rotation, which occurs when radio waves propagate through a magnetized plasma. 
Such a transformation may be written as r^) ((/>). Equivalently, the transformation of 
the Stokes parameters is R{)(2</>). A similar transformation generates a phase delay 
between the two components of E, and may be written as the rotation Yq{tp) . In 
optics a physical device that induces such a phase delay is called a compensator. This 
transformation arises in electronics when two signals traverse different cable lengths. 
Both r{)((/)) and r^('0) are unitary, and preserve the degree of polarization |S|/5o and 
the invariant interval 5inv 

Differential amplification or attenuation of the components of E provide examples 
of a non-unitary transformation. Consider the transformation 

g=f t ° ) (14) 

We may write this as g = ^QaQb bg(/?), where /3 = \\i{g}j/ga). The Stokes parameters 
transform as gagbBq{2(3). Note that this transformation does not preserve |S|/5o, and 
preserves 5inv only up to the factor gagt- 
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Next consider two orthogonal elliptically polarized waves with axial ratio tan^ and 
orientation 6 (Chandrasekhar 1960). 

/ cos cos y — i sin sin Y \ / sin^cos y — zcos^sinv \ ,,,, 

\ — sm a cos x ~ ^ cos y sm x J \ cos y cos x + ^ sm sm x / 

We may perform a change of basis by forming a matrix s with rows e\ and ej^. This 
matrix may be factored as s{6,x) = ru(x)r{)(^)- A transformation from a linear to a 
circular basis occurs for x = '?i"/4, and in the case where 9 = ir/A results simply in a 
cyclic permutation of the indices {q,u,v) — > {v,q,u) in equation |2|. The transformation 
s{9, x) rnay also be regarded as representing a receiver with two receptors sensitive 
to orthogonal forms of elliptical radiation. The process of reception constitutes a 
projection of E onto the receptors of the receiver, which is represented by matrix 
multiplication of E by s{9, x)- Clearly the choice of a rotation about the u axis followed 
by one about the v axis is not unique. One alternative specification of elliptically 
polarized radiation is given in terms of an orientation and phase delay as r{j{9)r^{(p). 
In optics this transformation is accomplished through a device known as a Babinet 
compensator (Born &: Wolf 1980). 

Now let us consider a receiver sensitive to two forms of elliptical polarization that 
are not necessarily orthogonal. This situation may arise in practice from imperfections 
in the construction of a receiver (Conway & Kronberg 1969, Stinebring et al. 1984). 
The transformation may be written as 



cos 9a cos Xa + i sin 9a sin Xa - sin 9a cos Xa + i cos 9a sin Xa 
sin 9b cos Xb + i cos 9i, sin Xb cos 9^ cos Xb — i sin 9f, sin Xb 



(16) 



where the two probes of the receiver are sensitive to elliptical radiation with axial ratios 
Xa, Xb and orientations 9a, 9b- For the case 9a = 9b and Xa = Xb equation |l^ simplifies 
to the unitary transformation s(^,x)- In the general case such a transformation does 
not conserve energy. With the definitions 



(^e = Oa+ Oh = Xa + Xb 

5e = 0a- Ob 5x = Xa- Xb 



(17) 



we may write c = c'r{)(cj0/2), where 



c' 



( cos{a^/2 + 5y,/2) cos 80/2+ - cos{a^/2 + 5^/2)sm5e/2+ \ 

i sin(cj;^/2 + 5^/2) sin 5e/2 i sm[a^/2 + 5^/2) cos 5q/2 

— cos{a-^/2 — 6^/2) sin Sg/2+ cos(cJ;^/2 — 6^/2) cos Sg/2+ 

y ism{a^/2 — 6^/2) cos 60/2 ism{a^/2 — 6^/2) sm6g/2 J 



(18) 



For a nearly orthogonal receiver with receptors sensitive to linearly polarized 
radiation, this matrix may be written to first order in these parameters as 
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c' = h'^l\50/2)h^y^ {5'x^/2)Y'^l\ay(^/2). Here the superscript (1) indicates that these 
transformations are first order in their arguments. 

These examples may be combined to form a model for the propagation and 
reception of radio waves. 

t = ,/^i^i^{P)v^{^i)hf {5e/2)hf {5j2)vf {a^/2)vi,{^^ 

(19) 

where $iono and $ism are the angles arising from Faraday rotation in the ionosphere 
and the interstellar medium, respectively, C is the angle between the frame of the 
receiver and that of the sky, and is the instrumental phase delay arising from 
differing electronic pathlengths. The equivalent transformation law for the Stokes 



parameters is simply obtained from equation 19 



S' = gagbBq{2f5)B.^{2^i)Bf{5e)^^^\5^)Yi^}\a^)^i,{^^^ S 

(20) 

The analysis of this model is simplified through the commutation relations in equation 
|l^. Amplifications and phase delays in the downconversion chain represent boosts and 
rotations with respect to the q axis, so it is easy to see that the order of amplifiers and 
relative electronics delays does not matter. Terms from individual components may 
simply be collected into the overall parameters (3 and Similarly, rotations about 
the same axis commute, and rotations about the v axis by the angles C, ^iono, and 
^ISM all have the same signature. Note that each of the parameters in this model may 
have its own time and frequency dependence. For example, $iono fluctuates in time 
as the ionospheric column density changes, and scales as from the cold plasma 
dispersion relation. Finally, calibration of a polarimetric observation is accomplished 



through the inversion of equations |19| or Naturally such an inversion requires a 



knowledge of the parameters in this model. 



5. Discussion 

The representation of polarimetric transformations in terms of the Lorentz group 
provide a simple context in which to analyze polarized radiation. This formalism is 
particularly relevant for the calibration of polarimetric data, and greatly simplifies the 
discussion from a qualitative standpoint. Propagation or instrumental effects that give 
rise to the rotations of equation ^ change the polarized component of the radiation S, 
thus obscuring the properties of the true, emitted light. Those effects that take the 
form of a boost transformation mix the total intensity So and the polarized component 
of the radiation S. Such transformations can have a particularly detrimental effect on 
polarimetric observations. In many astrophysical applications So is much larger than 
I S I , so that even boosts nearly equivalent to the identity matrix may completely corrupt 
the polarized flux. Observations that aim to detect very small polarized fractions, 
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such as the polarized component of the Cosmic Microwave Background radiation or 

the circularly polarized radio emission from Active Galactic Nuclei, are particularly 
vulnerable. For applications that require extremely high precision, the mixing of S 
into So can corrupt an observation. One such example has been seen in high-precision 
pulsar timing, where differential amplification and receiver cross-talk induce time 
dependent mixing of the pulse profiles, thereby modifying the total intensity profile and 
systematically shifting the times of arrival (Britton ct al., in preparation). For pulsar 
observations the invariant 5inv proves particularly useful, since an invariant pulse 
profile may be formed that is independent of propagation and instrumental effects. 
This invariant profile may then be used for timing or to investigate pulse variability 
that may be intrinsic to the pulsar. 

I thank Rene Grognard, Richard Manchester, Geoffrey Opat, and Matthew Bailes 
for useful conversations. I thank the referee for pointing out recent literature on the 
application of the Lorentz group to optics. 
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